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The two purposes of the paper are (1) to present a regularization of the self- 
^ ' field of point-like particles, based on Hadamard's concept of "partie finie", 

that permits in principle to maintain the Lorentz covariance of a relativistic 
field theory, (2) to use this regularization for defining a model of stress-energy 
^ . tensor that describes point-particles in post-Newtonian expansions (e.g. 3PN) 

of general relativity. We consider specifically the case of a system of two 
point-particles. We first perform a Lorentz transformation of the system's 
variables which carries one of the particles to its rest frame, next implement 
the Hadamard regularization within that frame, and finally come back to the 
original variables with the help of the inverse Lorentz transformation. The 
Lorentzian regularization is defined in this way up to any order in the relativis- 
tic parameter 1/c 2 . Following a previous work of ours, we then construct the 
delta-pseudo- functions associated with this regularization. Using an action 
principle, we derive the stress-energy tensor, made of delta-pseudo-functions, 
of point-like particles. The equations of motion take the same form as the 
geodesic equations of test particles on a fixed background, but the role of the 
background is now played by the regularized metric. 
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I. INTRODUCTION 



In recent years, the problem of the dynamics of gravitationally interacting compact 
objects in general relativity has received a lot of attention. This is due in part to the 
interest of the theoretical problem in its own, and in part to the ongoing development 
of laser-interferometric detectors for observing gravitational radiation. In the absence of an 
exact solution of the problem, one has recourse to successive post-Newtonian approximations 
(formal expansions in powers of 1/c). Within such approximations, it makes sense to model 
the compact objects with some "point-like particles", exactly as we do in a standard way 
within the Newtonian theory. However, the self-field of point-particles is infinite at the very 
location of a particle, and thus must be somehow regularized. The regularization is quite 
straightforward in the Newtonian theory, but it becomes non-trivial when going to high post- 
Newtonian approximations. Dealing with this problem, the present authors |1J developed a 
method for regularizing the infinite self-field of point-particles, which is based on the concept 
of "partie finie" , in the sense of Hadamard [^,[3] , of a singular function at the place of one 
of its singular points (see e.g. @-|7| for entries to the mathematical literature). We know 
that the Hadamard regularization yields the correct result for the equations of motion of 
two particles up to the so-called second and half post-Newtonian (2.5PN) approximation, 
corresponding to the order 1/c 5 beyond the Newtonian acceleration. Indeed, the problem 
has been completely solved at that order |@-@; notably some derivations make use of 



this regularization (e.g. fl2]J19|| ). In the present state of the art, we are concerned with 
the 3PN (or 1/c 6 ) approximation [EH-p^j. In fact, starting at this high post-Newtonian 



order, the regularization may become physically incomplete because of the appearance of an 
undetermined coefficient in the equations of motion ppj-pB]. 

The Hadamard regularization, investigated in is performed in a three-dimensional 
Euclidean space with generic point x e M. 3 , which is viewed as the spatial hypersurface 
labelled by t =const in a global coordinate system {t, x} covering the whole space-time. In 
particular, the regularization involves a spatial average, performed at t =const, over the 
direction of approach to the singularity. As such a regularization makes use of a preferred 
spatial hypersurface t =const, it is clearly incompatible with the framework of special or 
general relativity, which embodies a global Lorentz (or Poincare) frame invariance. Notably, 
we expect that the post-Newtonian equations of motion of point-like particles in harmonic 
coordinates (which we recall preserve the global Lorentz invariance) should exhibit at some 
stage a violation of the Lorentz invariance due to the latter regularization. The fact is 
that the breakdown of the Lorentz invariance due to the regularization occurs only at the 
very high 3PN approximation. Untill the 2.5PN order, it is sufficient to regularize within a 
preferred slice t =const of the harmonic coordinate system to obtain some Lorentz-invariant 



equations of motion [IS 



The first purpose of this paper is to define a regularization a la Hadamard that is 
compatible with the Lorentz structure of a relativistic field theory. This completes the defi- 
nition, proposed in ffl , of a specific version of the Hadamard regularization (based notably 
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on a particular class of pseudo-functions). To achieve this purpose, we shall simply perform 
the standard Hadamard regularization within the hypersurface that is geometrically orthog- 
onal, in the sense of the Minkowski metric, to the four-velocity of the particle. In separate 
papers [23L2l)|, we apply the latter "Lorentzian" regularization (together with the distribu- 
tional derivatives introduced in PJ), to the computation of the binary equations of motion at 
the 3PN order in harmonic coordinates, and find that, indeed, it permits the preserving of 
their Lorentz invariance (in some case at the price of adjusting some parameter). A different 
approach to the problem of incorporating the Lorentz invariance in the 3PN equations of 
motion consists of deriving a generic regularized dynamics, within the ADM-Hamiltonian 
formalism of general relativity, involving an arbitrary regularization parameter, and to de- 
termine this parameter uniquely by requiring the Lorentz invariance ||23|| . (See Section 2 in 
| 25| for a discussion on our point-mass regularization and its relation to p^ .) 

All-over the paper, we assume the existence of a preferred Minkowski metric, as selected 
for instance by the condition of harmonic coordinates in general relativity, with respect to 
which the trajectories of particles are represented by accelerated world lines like in special 
relativity. Most of our investigation is valid not only in the case of the gravitational field but 
also for any Lorentz-tensor field propagating on the Minkowski background. Furthermore, 
we shall define the Lorentzian regularization in a sense of formal expansion series in 1/c 2 ; 
so that, all the formulas in the paper will be given by some infinite series of relativistic 
corrections when c tends toward infinity. This is all what we need for the derivation of the 
equations of motion to the 3PN order [|24|,|25|. 

Since we are interested in the application to the motion of two particles, we shall define 
the regularization around one of the particles (say particle 1), and shall consider that its 
acceleration is purely due to particle 2. [However our definitions could be generalized to a 
system of N particles.] Notice that the particle 2 enters this regularization scheme through 
the Lorentz transformation of its own variables to the rest frame of particle 1, and the 
replacement of the acceleration of 1 in terms of the equations of the binary motion. In 
general, working at some given relativistic order, we shall need to know the equations of 
motion up to a lower order only, therefore giving us the possiblity of an iterative process. In 
this paper, we always assume that we know the relevant equations of motion at this order, 
and that these are Lorentz-invariant. 

Our second purpose is to derive an expression, compatible with the latter regularization, 
for the stress-energy tensor of point-like particles in post-Newtonian expansions of general 
relativity. Thanks to this regularization, we are able to give a sense to the value of the metric 
coefficients at the very location of the particle. Our basic assumption is that the matter 
action is the same as for test particles moving on a prescribed background gravitational 
field, except that the metric at the location of the particles is replaced by its regularized 
value in the sense of the (Lorentzian) regularization. From this assumption, we prove that 
the Dirac measure in the stress-energy tensor must be replaced by a certain generalized 
function defined by means of the Hadamard prescription. In the case of two particles (the 
generalization to N particles is immediate), we obtain 
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T P- lcle = , \ \ a Pf ^f=f + 1 - 2 , (1.1) 

where m% is the mass of the particle 1, and v f = (c, Vi) its coordinate velocity, i.e. vx = 
dyi/dt, y 1 = yi(t) being the trajectory parametrized by the coordinate time t (the symbol 
1^2 denotes the same expression but corresponding to the second particle). The notation 
[g po -} 1 means that the metric g p(J (t,x) is to be computed at the point x = yi(t) following 
the regularization (of course [g pa \i depends on the positions and velocities of both particles 
1 and 2). Note that the first factor in ( |1.1| ) is a mere function of time t. The second factor 
Pf ^ A ^^ yi ^ is made of a special type of partie finie delta-pseudo-function associated with 
the regularization (following the definition given in |[T|). It involves (minus) the determinant 
of the metric g pa , namely g, evaluated at the point (t, x), and a generalization Pf A(x — y x ) of 
the Dirac function defined in such a way that its action on a singular function yields the value 
of the function at the singular point in the sense of the regularization. Among the rules for 
handling the delta-pseudo-functions, we are allowed to write Pf ^ A( -*j7^ j = -^PfA(x — yi), 

whereas it is strictly forbidden to replace the latter quantity by [^^Pf A(x — yi). 

The stress-energy tensor ( |1 . 1| ) takes the same form as the one of test particles moving in 
a fixed background, but with the role of the background played by the regularized metric 
generated by the bodies. In particular, the equations of motion obtained from the covariant 
conservation of that tensor (V u T^ rticlc = 0), take the same form as the "geodesic equations", 
when considered with respect to the regularized metric. Our definition of the stress-energy 
tensor (|1 . 1|) constitutes a proposal, that we have found to be the most natural in the problem 
of the equations of binary motion at the 3PN order [0,0 , but that we have not proved to 
be generally valid to higher post-Newtonian orders (nor of course when considered outside 
a framework of post-Newtonian expansions). The tensor ( |1 . 1| ) appears to be a good candi- 
date for the characterization of point-like particles in post-Newtonian expansions of general 
relativity. 

The plan of this paper is the following. In Section 0, we recall from Q the material 
needed in the subsequent parts concerning the Hadamard regularization and the associated 
pseudo-functions. In Section [III], we investigate the formulas, needed to regularize, for the 
Lorentz transformation of some field point as well as two source points, and we define the new 
regularization around one of the particles as taking place within the instantaneous spatial 
hypersurface of the particle. In Section [TV], we give the formulas for this regularization at the 
level of the first relativistic correction 1/c 2 . Finally, in Section 0, we derive from an action 
principle our model of stress-energy tensor of point-like particles; the covariant conservation 
of this tensor leads to the equations of motion. 



II. HADAMARD REGULARIZATION 

To make the present paper self-contained, we shall review in this section the classic 
notions of the Hadamard regularization of singular functions and divergent integrals , 
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as well as the construction, by Blanchet and Faye |IJ , of a set of pseudo-functions associated 
with it. We follow closely the investigation of our previous paper |]J and employ most of its 
notation. A coordinate system {t, x} being given on space-time (for instance the harmonic 
coordinates used in Section IV), we consider some functions -F(x) defined on the spatial 
slice t =const, where x G R 3 denotes the position in the slice. We say that the function 
-F(x) belongs to the class T if and only if F is a smooth function on R 3 except at two 
isolated points yi and y2, and admits around each of these points the following power-like 
singular expansions. Denoting by r\ = |x — y x | the spatial distance to the point 1, and by 
n x = (x — yi)/r! the spatial direction of approach to 1, we assume that, for any N G N, 

F{*)= E r i AK) + °( r i) ■ (2-1) 

ao<a<N 

The coefficients if a of the various powers of r\ are smooth functions of the unit vector n.i, 
and the remainder tends to zero strictly more rapidly than r± when r\ — > 0. The powers 
a of ri in that expansion are assumed to be real, a G R, to range in discrete steps, i.e. 
a G (aj)ieN, and to be bounded from below, i.e. a < a for some ao G R. Similarly, we 
assume the same type of expansion around the point 2, 

ViVGN, F(x)= J2 r b 2 f b (n 2 )+o(r»), (2.2) 

b <b<N 

where r 2 = |x — y 2 | and n 2 = (x — y 2 )/T2- Thus, to each function F in the class T 
are associated two discrete families of indices a and b, and two corresponding families of 
coefficients i/ a (ni) and 2 /b(n 2 ), all of them depending on F. We shall refer to the coefficients 
if a for which a < (and similarly to 2 /b when b < 0) as the singular coefficients of F in the 
expansion when r\ — > 0. Since a > ao(F) and b > bo(F), the number of singular coefficients 
of F is always finite. 

The so-called "partie finie" in the sense of Hadamard 0|| of the singular function F at 
the location of the singular point yi is equal to the angular average, say i/ , of the zeroth- 
order coefficient, i/b( n i), in the expansion of the function when r\ — > we assumed in ( p.l|) ; 
namely 

= /» s / 17 /» (ni) • (2 - 3) 

where dQi = dQ(ni) denotes the solid angle element of origin y x and direction ni; the 
latter angular integration is performed within the coordinate hypersurface t =const. A 
crucial property of the Hadamard partie finie is its "non-distributivity" with respect to the 
multiplication, in the sense that 

(FG) l ^(F) 1 (G) 1 (2.4) 

in general. When applied to the gradient diF of a function F G J 7 , the definition ( |2.3| ) 
yields a useful formula which permits one to compute rapidly the partie finie of complicated 
expressions involving gradients: 



5 



(ftF)x = 3(^F) i . (2.5) 

Closely related to the concept of partie finie of a singular function is the definition of 
the partie finie (Pf) of the divergent integral J <i 3 x F. All-over this paper, we assume 
that the functions decrease fast enough at infinity (when |x| — > +00) so that the possible 
divergencies of integrals come only from the bounds located at the two singular points 1 and 
2. The "partie-finie integral" reads 0|| as 

Pf Sl , S2 I d 3 xF= lim If d 3 x F 

The integral in the right side extends over R 3 deprived from two closed spherical balls Bi(s) 
and B 2 (s) of radius s centered on the two singularities [thus Bi(s) and B 2 (s) are defined by 
ri < s and r 2 < s\. The other terms, which are defined by means of the partie finie in the 
sense of (|2.3|) , are chosen in such a way that the limit s — > exists. The notation 1^2 
indicates the same terms as the two previous ones but corresponding to the other point. The 
summation index a satisfies a < a < —3 (in particular the sum is always finite). Notice 
the two arbitrary constants si and s 2 which are introduced in order to adimensionalize 
the arguments of the logarithms in ( |2.6| ); the partie finie owns an ambiguity through these 
constants (hence the notation Pf sl)S2 ). The close connection between the partie finie of a 
singular function ( |2.3j ) and that of a divergent integral ( |2.(j| ) is most easily seen from the 
fact that 

Pf/ W ^K^) 1+ 1«2. (2.7) 

Unlike in the case of continuous functions, the (partie-finie) integral of a gradient is non-zero 
in general, and equal to the sum of the parties finies, in the sense of ( |2.3| ), of the surface 
integrals surrounding the singularities, in the limit where the surface areas tend to zero. 
This fact motivated the introduction and study in of a new derivative operator acting 
on JF, satisfying a property of "integration by parts" implying that the integral of any 
gradient is always zero. This operator generalizes for the class of functions JF the standard 
distributional derivative of Schwartz ||. 

Let us associate to any F G T a pseudo-function denoted PfF and defined to be the 
following linear form acting on the class T : 

WG E T , < PfF, G >= Pf / d 3 x FG , (2.8) 



where the right side is a partie-finie integral in the sense of ( |2.6| ); we use a duality bracket to 
denote the result of the action of the pseudo-function PfF on G. A fundamental definition 
adopted in 111], and motivated by the application to Physics, concerns the product of two 
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pseudo-functions, or of a function and a pseudo-function, which is the " ordinartf pointwise 
product in the sense that 



PfF . PfG = F . Pf G = G . PfF = Pf (FG) . (2.9) 

Thus, for instance, 

< PfF . PfG, H >= Pf J d 3 x FGH . (2.10) 

The product (|2.9|) chosen in [Q dictates most of the subsequent properties of the pseudo- 
functions, as well as their generalized distributional derivatives. (Refer to [[2~7| -|2~9|| for math- 
ematical treatises on generalized functions and distributions.) In particular, the derivatives 
do not in general satisfy the Leibniz rule for the derivation of the product, although they 
satisfy it in an "integrated sense" , according to the rule of integration by parts. 

The Riesz delta-function, given for e > by e 5(x) = |x| e ~ 3 , tends, in the usual 
sense of distribution theory, towards the Dirac measure when e — > 0. When considered with 
respect to the singular point yi, the Riesz delta-function allows us to define a useful element 
of our class, 

A(x)^(x- yi ) = ^rr 3 e T. (2.11) 

Therefore it is possible to associate to e 8\ (for any e > 0) the pseudo-function Pf £ Si following 
the prescription (|2.8|) . Applying the limit e — > 0, we obtain 

lim < Pf £ 5i, F >= lim Pf / d 3 x ^F = (F), , (2.12) 

e^O e^O J 



where the value of F at the point 1 in the right side is defined by the prescription ( |2.3|) . This 
motivates us for introducing a new pseudo-function, we shall call the delta-pseudo-function 
Pf#i, as the formal limit of the pseudo-functions Pf £ 5i when e — > 0. By definition, 

VFG J, < Pf5i, F >= (F) 1 . (2.13) 

Clearly, the delta-pseudo-function generalizes the notion of Dirac distribution 8\ = 
5(x — yi) to the case where the "test" functions are singular and belong to the class T . 
Extending the definition of the product Q2.9|) to include the delta-pseudo-function we pose 

Pf F . Pf 6i = F . Pf 8 X = Pf (FSt) , (2.14) 

as well as, for instance, 

Pf . PfG = Pf {FSx) . G = Pf(FG5i) . (2.15) 
The new object Pf (F5\) in ( p . 1 4| ) - ( pTT5| ) has no equivalent in distribution theory; it satisfies 



VGeF , < Pf (F5 1 ), G >= (FG) 1 . (2.16) 
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We notice for future reference that a consequence of the "non-distributivity" of the 
Hadamard partie finie [see ( j2.4| )1 is that 



Pf(F5x) + {F\ Pf5 1 . 



(2.17) 



We are not allowed to replace a singular function that appears in factor of the delta-pseudo- 
function at the point 1 by its regularized value at that point. 

The derivative of the delta-pseudo-function Pf^ was constructed in [[j]. As it turns out, 
it takes the form of an "ordinary" derivative : <9j(Pf5i) = Pf(<9j5i); due to the presence of 
the delta-pseudo-function, there are no distributional terms associated with it. We have 
simply (from the rule of integration by parts), 

VF G T , < ^(PWi), F >= — < PfJx, d,F >= -(d i F) 1 . (2.18) 

The differentiation of the more complicated object Pf(F5i) proceeds in the same way: 

VGgJ, < di\Pf(F5i)], G>=- < Pf (FSt), diG >= -{Fd i G) 1 . (2.19) 

we can write for the intrinsic form of this 



Note that, as a consequence of the identity 
object: 

9i[Pf (F<fi)] = Pf 



m 3 M 



(2.20) 



Because the derivative of the delta-pseudo-function is equal to the ordinary one, the Leibniz 
rule for the derivative of a product happens to still hold. For instance, in the case of the 
product of Pf(F5i) with some pseudo-function PfG, we have 



<9;[Pf(F5i). PfG] = di[Pf(F5i)]. PfG + Pf ^(PfG) . 
The proof uses the combination of (|2.15|) and ( |2.19|) . 



(2.21) 



III. LORENTZIAN REGULARIZATION 

To define a Lorentzian regularization a la Hadamard (based on the investigation of 
and on Section [TT]), we now need to specify in a precise way the dependence of a function 
F(x) in the class T on the "source" variables at the coordinate time t of a global frame 
{x, t}. We assume (as everywhere else in this paper) that we are working at some given finite 
order in a relativistic or post-Newtonian approximation. Up to a given order, we can choose 
as the source variables the two trajectories yi(i) and y2(t) in the frame {x, t}, and the two 
coordinate velocities Vi(i) = dyi/dt and v 2 (i) = dy 2 /c?t (the trajectories of the particles are 
time-like world lines in Minkowski space-time). Indeed, it is legitimate to assume only the 
latter source variables because, up to a given post-Newtonian order, we can order-reduce the 
accelerations and all derivatives of accelerations by means of the equations of motion of the 
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particles up to the appropriate accuracy (in general the precision of the equations of motion 
needed for this order-reduction is one order less than the given post-Newtonian order at 
which we are performing a calculation). Of course, we are assuming that these equations of 
motion are known (they are known presently to the 2.5PN order |TB|JT^JTP[] , and the general 



motivation of this work is to get them up to the 3PN order [2"2",pi,2l)|. Thus, we assume 
that the function F G T really reads 

F(x,t) = F[x;y 1 (t),y 2 (t);v 1 (t),v 2 (t)] . (3.1) 

We denote with the same letter F, by a slight abuse of notation, the function of the field 
point (x, t) and the functional of the field point and source variables in the right-hand-side. 
For definiteness, we assume that the two trajectories are smooth functions of time, i.e. yi, 
y 2 G C°°(]R 3 ), and that F is a smooth functional of the two velocities v 1; v 2 (see also Section 
IX of |I| for details about our assumptions). By (|3.1|), we mean that the dependence of F on 
the coordinate time t is through (and only through) the two instantaneous trajectories yi, 
y 2 and velocities v 1; v 2 . Note also that it is implicitely assumed with our notation (|3.1|) that 
the function F depends locally on time t (no dependence over the trajectories and velocities 
at some time earlier than t for instance). Furthermore, very often in applications, we shall 
find that the dependence of F on the spatial position x appears only via the two spatial 
distances to the source points, ri(t) = x — yi(t) and r 2 (t) = x — y 2 (t). In this paper, we 
shall generally suppose, in order to simplify the presentation, that this is the case; namely, 
the function F, as a functional of the source variables, is 

F(x,t) = J P[r 1 (t),r a (t);v 1 (t),v 2 (t)] . (3.2) 

The hypothesis ( |3.2|) does not constitute a very severe restriction. The extension to the 
more general case (|3.1| ) is generally straightforward; moreover, (|3.2|) is always verified in 
the problem of the post-Newtonian equations of motion of binary systems. In this section, 
we shall define the Lorentzian regularized value of the function F at the location of the 
singularity 1, by contrast to the non- invariant regularized value defined by ( |2.3| ) within 
the "global" coordinate hypersurface t = const. We shall denote by [F] 1 the new Lorentzian 
regularization of F at the point 1, defined within the instantaneous rest frame of the particle 
1 at t' =const [in constrast with the notation (F) 1 used in ( |2.3| ) for the old regularization]. 
In addition, we shall introduce a delta-pseudo-function denoted by Pf Ai associated with the 
new regularization [similarly to the delta-pseudo-function Pf^ which was defined in ( p. 13 ) 
in the case of the old regularization]. 



A. Lorentz transformation of the source variables 

In this paper, it is sufficient to consider only those homogeneous proper Lorentz transfor- 
mations which change the velocity of a global inertial frame {x^} = {ct, x}. More specifically, 
let us consider the Lorentz boost 
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a/" = A"„(V) re" , (3.3) 
where the Lorentz matrix A^(V), depending on the constant boost velocity V, is given by 

A°o(V)= 7 , (3.4a) 
A 4 (V) = -7^ , (3.4b) 

A° J .(V) = -7^, (3.4c) 
J c 

A W = ^ + 7TiV- (3 - 4d) 

We indifferently denote the components of the boost vector by V = (V 1 ) = (Vi) (spatial 
indices i,j = 1, 2, 3). The Lorentz factor 7 reads 

7 = -p=L= , (3.5) 

with V 2 = 5ijV l Vi (of course |V| < c). The inverse transformation is x u = A^(V)x ,fJi where 
the components of A ^(V) = r]^ p r] ua A p a (y) are obtained directly from (|3.3|) by changing 
V — > —V. The choice of sign made in the Oi components of the boost (|3.4j) is such that a 
particle which has velocity V at time t in the frame {x^ 1 } is at rest in the frame at 
time t' . 

We introduce on one side the space-time event Q, which represents for us a "field" point 
located outside the two world lines of the particles, and on the other side the space-time 
events Pi, Mi and P2, M 2 , which are "source" points, lying respectively on the world lines 
of the particles 1 and 2 (see below for their definition). The coordinates of the event Q are 
(t, x) in the frame and (t', x') in the frame {x^}. Sorting out the spatial and temporal 
indices in fl3.3|), we have 

ct' = cA° t + A jX j , (3.6a) 



x 



11 = cA\t + A^x j . (3.6b) 



The points Pi and P2 are now defined as the two events that are located on the trajectories 
of the particles and are "simultaneous" with the event Q in the frame i.e. that 

belong to the same spatial slice t =const as Q. The coordinates of Pi and P 2 in {x 11 } are 
denoted by (£, yi) and (t,y2) respectively, the two trajectories yi = yi(t) and y 2 = y2(£) 
being parametrized by the coordinate time t in that frame. On the other hand, in the new 
frame {x /f1 }, the coordinates of Pi and P 2 are (r^z^) and (t2,z 2 ). Evidently, the primed 
coordinates are related to the unprimed ones by the Lorentz boost ( |3.3|) , so that 



ct[ = cA° t + A° jy { , (3.7a) 
^ = cAV + A^, (3.7b) 
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in the case of the event Pi [where y{ = y{(t), y 2 = yl(t)\, and 

cr^ = C AV + A ^, (3.8a) 
z' 2 l = cA\t + A l jy i, (3.8b) 

in the case of the event P2. In the new frame the source events that are simultaneous 

with Q are not Pi and P2, but some other events Mi and M2, whose coordinates in the 
primed frame are thus (t ; , y[) and (£', y 2 ); the coordinate time t' is the same as that of Q in 
the primed frame, and the spatial coordinates are the trajectories of the particles y[ = y[(t') 
and y' 2 = y' 2 (t f ) which are labelled by t' in the new frame. Let (ti,Zi) and (t 2 ,z 2 ) be the 
coordinates of Mi and M 2 in the original frame By definition, 

ct' = cA° ri + A° jZ i , (3.9a) 

y[ l = cA\T 1 + A^z{, (3.9b) 

ct' = cA° r 2 + A^zi , (3.9c) 

y' 2 l = cA\r 2 + A^zi, (3.9d) 

where y[ l = y'i(t') and y 2 = y' 2 (t'). Beware of our notation, where t[ (for instance) is the 
time coordinate of Pi in {x ,f1 } while r± is the time coordinate in {x^} of the different event 
M 1. Since the events Mi and M 2 are located on the world lines of the particles parametrized 
by yi{t) and y 2 (t) in it is clear that at time T\ in that frame their coordinates are 

related to the trajectories by 

z i = yi( r i) , (3.10a) 
z 2 = y 2 (r 2 ) . (3.10b) 

Similarly, from the fact that Pi and P 2 are also on the world lines, which write as y'i(t') 
and y' 2 (t') in the frame {x /f1 }, we deduce that their coordinates in {x'^} satisfy 

z i = yi( r 0> (3.11a) 
z' 2 = y 2 (r^) . (3.11b) 

By eliminating t' from the equations (|3.6a|) and ( |3.9a ) we immediately obtain 



cA\(r 1 -t)=A\(x i -z\), (3.12) 



or, equivalently, taking also into account (|3.4|) 



n-t = -^V.(x-z a ), (3.13) 



e 2 



where the usual Euclidean scalar product between (boldface) vectors is denoted by a dot. 
With the help of the latter formula for expressing n, we can re-state the belonging of zi to 
the particle world line at time t\ [see (|3.10a )1 as 
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zi = yiU-^V.(x- Zl ) ) . (3.14) 



c 2 



Recall that zi is the spatial coordinate in the old frame of the event Mi which is simultaneous 
with the field point Q in the new frame. Clearly, the equation ( P-14j) determines the vector 
zi as a function of the coordinates (t, x) of the field-point event Q (see the appendix). Here, 
let us view zi as a "vector" field zi(x), solution of ( P-14j ), lying in the three-dimensional 



space t =const. It is evident from Q3.14 ) that the function Zi(x) admits a fixed point at 
Yi = yi(*)> i n the sense that 

zi(yi)=yi- (3-15) 



Unless specified otherwise [like in ( p. 14 )1, the notation yi always means yi(i). The mathe- 



matical justification of (|3.15|) is the following. From the fact that the world line of the particle 
is time-like we can write, for any instants t and t, the inequality \yi(t) — yi{t)\ < c\t — t\. 
Hence, applying the definition ( p,14| ), we find that our function Zi(x) obeys, for any positions 
x and x, the further inequalities |z x (x) — Zi(x)| < i|V.(x— x)| < ^-|x — x|. Now recall that 
— < 1, so the latter inequalities mean exactly that the function x — ► zi(x) is a contracting 
application with respect to the usual Euclidean norm (i.e., it satisfies the property of Lips- 
chitz with a ratio k — — strictly less than one). Therefore, by the theorem of Picard, the 
function admits a unique fixed point, which of course is nothing but yi. (Besides, at the 
location of the fixed point, we have T\ = t.) 

In this paper, we establish the general solution of the equation ( |3.14|) in the form of 



an infinite (post-Newtonian) power series in 1/c 2 . We shall not discuss the convergence 
properties of this series and simply employ it to define the regularization up to any relativistic 
order. This is sufficient for the application to the problem of the equations of motion of 
particles in the post-Newtonian approximation. The general solution of ( 3.14j ), as determined 
in the appendix, reads 



+oo 



zi = yi 



n=l 



n—1 






"(V. ri r Vl " 



+ CV-rJ-v, , (3.16) 



with shorthand notations yi = yi(t), = x — yi(i) and vi = vi(£). The many derivatives 
d/dt in the right side are partial time derivatives with respect to the coordinate time t, 
the spatial coordinate x being held constant. They act on ri through the trajectory y^: 
we have dv\/dt = — vj or d(V .r\) /dt = — V.vi for instance. On the other side, they act 
of course on velocities and (derivatives of) accelerations: thus dvi/dt = ai, d&i/dt = b 1; 
dhi/dt = Ci, and so on, where a 1; bi, Ci represent the acceleration, and its first and second 
derivatives (in these cases the partial derivative is a total derivative, e.g. dwi/dt = a^). 
Thus, to high post-Newtonian order, ( |3.16| ) contains many accelerations and derivatives 



of accelerations, but it is understood that this formula is order-reduced, consistently with 
the post-Newtonian order; i.e. all accelerations and derivatives of accelerations are to be 
replaced by the functionals of the positions and velocities deduced from the equations of 
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motion. Combining ( |3.13|) and ( |3.16|) , we easily find the corresponding solution for the time 
coordinate Ti, 



+00 

n=l 



c 2n n\ 



n-l 



(v.n)» 



(3.17) 



[Of course, since V is a constant vector, it could be as well put outside the partial time 
derivative operators in both (|3.16| ) and ( |3.17| ).1 Finally, equations ( |3.16| ) and ( |3.17| ) determine 
completely the space-time event Mi. From them, we can recover directly the fact that when 
x = yi (at the fixed point) then Zi = yi and T\ = t: there are in the right sides of 
both relations n — 1 partial time derivatives acting on a term that involves the nth power 
(V.ri) n , so that at least one of the scalar products V.ri is left un-differentiated, and makes 
the sums in (|3.16| )- (|3.17| ) vanish when ri = 0. Replacing both z x and T\ as given by the 



infinite post-Newtonian series back into the relation ( 3.10a ), expressing both sides of the 
resulting equation as the same type of post-Newtonian series with the help of a formal Taylor 
expansion when c — > 00, and finally equating all the coefficients of these two series, yields 
an interesting mathematical formula relating together some sums of products of derivatives. 
This formula is derived in the appendix (where we present also a direct proof of it). Notice 
that the same reasoning as before can be done on the coordinates of the event Pi in the 
new frame, that we find to be given by 



+00 ^ / a s n-l 

c 2n n\ \dt 



n=l 

+00 ^ ( Q \ n ~ ' 
Tl=t + Y1 (foi 

n=l x 



(v.ri)v; 

v.r;) n 



(3.18a) 
(3.18b) 



where y[ = y'i(f), r[ = x' — y'i(t') and v[ = v'i(t'). Evidently, the result Q3.18D can also be 
deduced directly from ( p.!6| )-( p.l7| ) by changing V into —V and replacing all the non-primed 
variables by the corresponding primed ones. 

We are now able to give all the transformation laws of field and source variables between 
the frames and {x ,f1 }. Of course, from ( |3.6| ), the transformation of the field variables 
is the standard Lorentz one, 



t' = 7 (t-i(V.x)) 



x - 7V [t 



7 



c 2 7 + 1 



(3.19a) 
(3.19b) 



Concerning the source variables, we are interested in the expressions of the new positions 
yi(i'), y^') an d velocities v[(t'), v' 2 {t') in the new frame at time t'. These are straightfor- 
wardly obtained from inserting the results ( |3.16| ) and ( p. 17] ) into the equations (pT9|), as well 
as the similar results corresponding to the point 2. We find, for trajectories, 
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y'i = yi - 7V ( t 

+00 

+ E 

n=l 



1 7 
c 2 7 + 1 



c 2n n! V dt 



y 2 = y 2 - 7V ft 

+00 / 

+ E 



1 7 

c 2 7 + 1 

d\ n— 1 



V.x 

V.r, 
V.x 



Vi 



7 



7 + 1 



V 



(3.20a) 



n=l 



c 2n n! \<9t 



(V.r a ) 



v 2 



7 



7 + 1 



V 



(3.20b) 



By subtracting the latter equations (|3.20|) to x' as given by ( 3.19b ) we obtain the spatial 
distances = x' — y[(t') and r' 2 = x' — y' 2 (t') as 



ri 



E 



' c 2n n\ \dt J 



n-l 



(V.n) n Vl 



7 



7 + 1 



v 



r 2 



-E 

n=l 



c 2n n\ \dt J 



n-l 



(V.r 2 ) n v 2 



7 



7 + 1 



V 



(3.21a) 
(3.21b) 



These relations will play the crucial role in the definition of our Lorentzian regularization. 
Of interest also is the expression of the relative distance between the two particles, i.e. 
y'i2 = y'i - y'2 = r 2 - r i S iven b Y 



yi2 = yi2 



E 

n=l 



c 2n n\ \ dt 







n-l 



(V. ri ) n Vl 



7 



7 + 1 



V - (V.r 2 )" v 2 



7 



7 + 1 



V 



(3.22) 



Finally, we compute the expressions of the coordinate velocities v^(t') = dy^/dt' and v' 2 (t') = 
dy 2 /dt' in the new frame. They follow immediately from the law of transformation of the 
time derivative : d' t = jd t + r )V l d i) and we obtain 



/ = i i^(-r / ay 

' n=l N 7 



4 = -v 2 - V 

7 



1 +00 

^E 



7 * — ' c 2n n\ V dt 

n=l 



<9 



(V.r x ) n ( v, 
(V.r 2 ) 



7 



v 2 - 



7 + 1 

7 

7 + 1 



V 



V 



(3.23a) 
(3.23b) 



Notice that although the velocities v^(t') and v 2 (t') are some mere functions of the coordinate 
time t' in the new frame, they depend, when expressed in terms of quantities belonging to 
the old frame, on both time and space coordinates t and x. This is obvious because by 
changing the space coordinate x of the field point Q while keeping t =const we change the 
time coordinate t' of the source events Mi and M 2 and therefore the values of their particle 
velocities (soon as the trajectories are accelerated). This fact is important and has to be 
taken correctly into account in the regularization process defined in the next subsection. 
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The inverse formulas are obtained in the same way by substituting ( |3.18| ) into the inverse 
of ( |3.7| ). They correspond of course to changing V into —V and replacing everywhere the 
un-primed labels by primed ones. We find, for the spatial distances and velocities, 



I"! 



r 2 



vi 



v 2 



V.r' 2 f 



r > _ V— f — 

r l c 2n n \ I fa, 

n=l x 

r 2 c 2n n l [Qfl 

n=l v 

' 71=1 X ' 



v. r ;r v; + 



7 



Vo + 



7 + 1 

7 

7 + 1 



V 



V 



(V.r'O* vi + 



(v.r' 2 ; 



7 



7 + 1 

7 

7 + 1 



V 



V 



(3.24a) 
(3.24b) 
(3.24c) 
(3.24d) 



B. Definition of the regularization 

Let us consider a function F belonging to the class T and being at the same time a 
scalar under Lorentz transformations, i.e. -F(x, t) = F'(x.',t'). More precisely, we restrict 
ourselves to the case of a dependence on x only via the distances 17 and r 2 [cf ( |3.2|) 1; this 
means 

F[ ri {t)Mt)-Mt)Mt)] ^Fy^y^yy^y^v] , (3.25) 

where we use the same slighly abusive notation as in ( p.2| ), with addition, in the right side, 
of the explicit mention of the dependence over the boost vector V. All the variables in both 
frames {x 11 } and {x ;/1 } are related to each other by the formulas developed in the previous 
subsection. The regularization process goes as follows. 

(I) Starting from F[r 1; r 2 ; v 1; v 2 ] defined in the frame {x^}, we first determine the new 
functional F'fr^, r' 2 ; v' 1; v' 2 ; V] in the boosted frame {x'^}. To do so, we replace all the 
variables 17, r 2 , Vi, v 2 by their expressions in terms of the new ones r^, r' 2 , v^, v' 2 as given by 
the formulas ( |3.24j ), in which it is understood that all the accelerations are order-reduced up 
to some given specified post-Newtonian order. Performing all the necessary post-Newtonian 
re-expansions to that order, we indeed obtain in that way (since F is a Lorentz scalar) the 
new functional F' of the new distances r' 1( r' 2 and velocities v' 1; v' 2 . In addition, F' depends 
as expected on the constant V which is yet un-specified at this stage. 

(II) We compute the Hadamard regularization of F' at the point 1 following exactly the same 
rules as defined in (|2.3j ), but in the boosted frame {x'^} (in particular, within the coordinate 
slice t' =const). In words, we perform the expansion of F' when the spatial distance r[ 
tends to zero, and obtain the same type of power-law expansion as in (|2.1| ) [since the form 
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of the relations Q3.24Q shows that the structure of the expansions in both frames must be 
the same]. However, we get some primed functional coefficients \f' a that differ from the 
un-primed coefficients \f a appearing in ( |2.1| ). The boost vector V is simply held constant 
in the process. Thus, WN e N, 

F'[r;,r' 2 ;v;,v^V]= £ r'l ft (n' i; y' 12 ; v[, v 2 ; v) + o (r'f ) , (3.26) 

a <a<N 

with the notation r'\ = |x' — y\\, n[ = (x' — y'\)/r'\ and y' 12 = y[ — y 2 . (The fact that the 
coefficients \f' a depend on y' 12 instead of the two individual trajectories y[, y' 2 is due to our 
restriction that F' depends on x' via the distances r[ and r' 2 ; also, the accelerations depend 
on the relative distance y' 12 .) Now, like in ( |2.3| ), we pick up the zeroth-order coefficient in 



the r^-expansion ( |3.26| ) and average over the angles. This defines a certain functional of the 



separation vector y' 12 , the velocities v' 1; v' 2 and the boost velocity V, 

/ofy^'i'^ V ) = /^T {o^jy^v^v) . (3.27) 

We insist that the angular average is performed in the new frame, within the spatial hyper- 
surface t' =const; in particular, the solid angle element in (p.27 ) is the one associated with 



the unit direction in that hypersurface: dQ\ = dQ(n'i). Here again, V is considered as 
a simple constant "spectator" vector during the average. 

(Ill) We impose that the new frame is actually the rest frame of the particle 1 at the event 
Pi. Recalling that the Lorentz boost ( |3.4|) brings a particle with velocity V in the frame 
at rest in the frame {x" 1 }, we see that we must choose 

V = vi(t) . (3.28) 

We come back to the original variables in the un-primed frame by using the transformation 
laws ( |3.22| )-( |3~23| ), in the limit where the field point x tends to the source point yi(t) (because 
we are located at the event Pi), with V = vi according to ( |3.28 ). Note that, in this limit 



r x — > 0, the coordinate time t' of the event Q in the primed frame is equal to the coordinate 
time t[ of the event Pi. It is important to realize that both the computation of the limit 
when ri — > and the replacement of the vector V by (|3.28|) are to be done after performing 



the many partial time differentiations in ( |3.22|) -( |3.23|) . Consider first the primed variable 
y' 12 , which is given by ( |3.22| ) where we apply the replacement ri = (as well as V = vi). 
In ( |3.22|) the n — 1 partial time derivatives acting on the term proportional to (V.ri)" - will 
clearly lead to zero in the limit ri = 0; indeed, by an argument met previously, there are 
not "enough" derivatives to make a non-zero contribution. So the variable to be used when 
coming back to the original frame is 



(3.29) 



r2=yi2 
V=vi 
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As indicated by the notation one must implement the replacements of r 2 by yi2 (this is 
equivalent to i"i = 0) and of V by vi after the n — 1 time differentiations. In the case of the 
primed velocity of the particle 2, given by (3.23b), we simply have 



7 



-V 2 -V 



7 *■ — ' c 2n n\ 

71=1 



(V.r 2 ) 



v 2 



7 



7 + 1 



V 



(3.30) 



r2=yi2 
V=vi 



The formulas ( |3.29| ) and (|3.30| ) define, after order-reduction of the accelerations, some func- 
tionals y' 12 [yi2; Vi, v 2 ] and v 2 [yi2; Vi, V2] that we use for coming back to the initial frame 
Clearly, the primed velocity v[ of the point 1, at which we perform the regularization, 
deserves a special treatment. From ( |3.23a| ) we obtain 




7 



7 + 1 



V 



(3.31) 



n=o 

V=vi 



Here, there are n time derivatives which is a priori enough to make a contribution. The 
only possibility is to differentiate successively each of the n factors V.r 1; yielding for each of 
the terms in the sum n\ identical contributions. Hence, we arrive at a much simpler series, 



7 



- Vl _ V 



+00 

' n=l 



V.vi 



Vl 




(3.32) 



V=vi 



which can now easily be summed up. The result is 

-vi - V + -L^V S 

/ 1 7+1 cr 



1 



V.vi 



(3.33) 



V=vi 



from which we immediately deduce that the primed velocity of the particle 1 must be zero, 

vi = . (3.34) 

This is of course the expected result because the boost velocity was chosen to be equal to 
the instantaneous velocity of the particle 1 in the un-primed frame at the instant t; however, 
the details of the above proof constitute a necessary consistency check of the formulas. 

(IV) The choice of boost vector V = v 1; together with the equivalent statement that v[ = 0, 
as well as the expressions (|3.29| ) and (|3.30|) defining the two functionals y' 12 [yi2; vi, v 2 ] and 
v 2 [yi2; vi, v 2 ], are put into ( |3.27|) , which gave the result i/q of the spherical average in the 
Hadamard regularization performed in the primed frame. Therefore, the regularized value 
of F at the point 1 is defined by 



l F ]i = /o(yi2[yi2;vi,v 2 ];0,v 2 [yi2;vi,v 2 ];vi 



(3.35) 
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The new regularization [F] 1 acts, like the old one (F) 1 , as a certain functional of the relative 
distance yi 2 and the velocities vi, v 2 . However, in generic cases, [F] 1 differs from (F) 1 by 
relativistic terms at least of the order 1/c 2 [we investigate in Section |IV| the exact relation 
between both regularizations to the first relativistic order 1/c 2 ]. In the problem of the 
post-Newtonian equations of motion, we have found P2"|j2"4[] that the new regularization [F] 1 
adds some extra terms to the acceleration computed using the regularization these 
new terms are of order 3PN and manage to make the 3PN equations of motion invariant 
with respect to Lorentz transformations. Indeed, with the regularization (F) 1 the Lorentz 
invariance of the equations of motion would be broken at the 3PN order. Finally, let us 
introduce as we did in Jl| (see also Section [H]) a delta-pseudo-function associated with 
the new regularization [F] 1 . By definition, the "Lorentzian" delta-pseudo-function denoted 
Pf Ai [to contrast with the non- invariant one PfSi defined by (|2.13|) 1 is such that 

VF G J 7 , < Pf Ai, F >= [F] x , (3.36) 

where the right side is given by the new regularization (|3.35| ). By definition, we have in 
the case of the new regularization the same laws for the multiplication as in Section ||, for 
instance 



FG . Pf A x = Pf(FAi) . G = Pf (FAi) . PfG = Pf (FGA X ) , (3.37) 

where the pseudo-function Pf(FAi) is defined by 

VGG J, < Pf (FAi), G >— [FG] 1 . (3.38) 

This pseudo-function Pf(FAi) is at the basis of our proposal for the stress-energy tensor 
of point-particles in Section |V|. And, like in the case of Pf(FSi), we are not allowed to 
replace this pseudo-function by the product of the regularized value of the function times 
the delta-pseudo-function, namely 

Pf(FAx) ^ [F], Pf Ai . (3.39) 

The derivatives of Pf Ax and Pf(FAi) are constructed in the same way as for the original 
regularization in Section |I[ Therefore, 

VG G T , < $[Pf(FAi)],G >= - < Pf(FAi), 9*6? >= -[Fd i G] 1 . (3.40) 



However, the identity ( |2.5| ) is not valid in the case of the new regularization, so we do not 
have a result similar to (|2.20|) [see ( [4.1 3D for the equivalent of (|2.5|) at the first relativistic 
order]. For the product of Pf (F5\) with some PfG, the Leibniz rule holds: 

di[Pf (FAi). PfG] = di[Pf(FAx)]. PfG + Pf(FAi). ^(PfG) . (3.41) 

This is a consequence of the definition (|3.40|) and the law ( |3.37|) . 
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IV. THE REGULARIZATION AT THE FIRST RELATIVISTIC ORDER 



At this point, it is instructive (and useful in practice) to present the complete formulas 
that define the Lorentzian regularization [F] 1 at the level of the first relativistic corrections 
1/c 2 , i.e. neglecting all the terms of order 0(l/c 4 ). [Notice that, consistently with Section fT!\, 
we must consider that the boost vector V itself is of order 0(1), so that for instance the factor 
V 2 /c 2 really represents a small relativistic correction of the order 0(l/c 2 ).] Furthermore, 
we shall obtain at this 1/c 2 level a formula linking the new regularization [F] 1 to the old 
one (F) 1 . Like in Section |HI| , we assume that the function F depends on x through the two 
distances ri(£) and r 2 (t) only; this implies a relation between the partial derivatives: 



0F + — + — 







(4.1) 



(where d. t = d/dx 1 ). We suppose also that F is a Lorentz scalar, cf ( |3.25| ). 

We follow the general specification for the regularization in Section |T|. We first express 
the vectorial distances r 1; r 2 and velocities v 1; v 2 in the boosted frame {x >f1 } using the 
transformation formulas ( |3.24| ) restricted to the order 1/c 2 . For the distances, we get 



ri = ri-i(V.ri)|Vi + iV 



r 2 



r 2 - V 



V.rV 



V2 + |V 



o 



o 



(4.2a) 
(4.2b) 



The relative distance yi 2 = r 2 — r x reads as 



yi2 = y i2 



1 



•Kvy^v + cv-riM- 



;V.r' 2 )v 



+ 



(4.3) 



while, for instance, the relative separation r 12 = |y 12 | is 



ri2 - r[ 2 1 + - -|(V.n' 12 ) 2 + i(V.n;)(v',n' 12 ) - -£(V.n 2 )(v 2 .n' 12 ) 



where = v' x /r' X) n' 2 = r' 2 /r 2 and n' 12 = y' 12 / r 'i2- F° r ^ ne ^ wo velocities, we find 



(4.4) 



Vl = vi + V + - ( [-|V 2 - V.vi] v; - KV.v'jV + (V.ri)ai) + 



v 2 = v' + V + 



-iV 2 -V.v 2 



v 2 -i(V.v 2 )V + (V.r 2 )a 2 )+0 



(4.5a) 
(4.5b) 



where the two accelerations a' x and a 2 are to be replaced, consistently with the approximation, 
by their Newtonian values : a; = -^n' 12 + (£) and a' 2 = ^n' 12 + (^). [Notice that 
in Section |TTT| the regularization has been defined regardless of the type of special-relativistic 
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interaction involved; in the case of electromagnetism, for instance, we should simply replace 
the accelerations by their Coulombian values in ( f4.5| )-l 

Next, we substitute the expressions ( [4.2| ) and ( fOS| ) into the scalar function F[r 1; r 2 ; v 1; v 2 ] 
and perform the expansion to the first order. The result is the scalar function 
F'fr^, r' 2 ; v^, v 2 ; V] in the new frame; thus 



F'[r[, r' 2 ; v' 1; v' 2 ; V] = F[r[, r' 2 ; vi + V, v' 2 + V] 



+^(V.r' 1 ) v'\ + \V 



M + ? (v - r2 



v'l + 



l dF 



+- 



+- 



-|v 2 -v. v ; 



iv 2 -v.v 2 



§(VV 1 )V 4 +(V.i' 1 )a'i 



^-I(V.v 2 )V" + (V.r' 2 K 



dF 
dv\ 
dF 



2, M + 



(4.6) 



where we have used |£ = — Jrr and 14- = —Jrr- Note that, to this order, the partial 
derivatives in (|4.6|) can be evaluated at the primed values r' 1; r' 2 and v[ + V, v 2 + V, or 
equivalently at the non-primed ones rj, r 2 and vi, v 2 . Now we pick up in the new frame 
the term of zeroth order in the expansion when r[ — > 0, and perform the angular average 
with respect to the direction n^. This yields the functional of the variables y' 12 , v[, v 2 and 
V which has been defined in ( |3.27p . Since these operations of expanding and averaging 
represent nothing but the Hadamard regularization in the old sense of ( |2.3| ), we can denote 
them by using the parenthesis appropriate for this regularization. Therefore, 



(o (y' 12 ; vi, v' 2 ; v) = ( F[ ri , ri + y' 12 ; v ; + V, v' 2 + V] 



+- 



+- 



+^(V.r 1; 



1 g + ^(v. ri + vy 12 



_I V 2 - V.v; v'[ - i(V.vi)V* + (V.r x )a'i 



A + \v l 

dF 

dv\ 



dF 



-v 2 -v.v; 



v n 2 -Uv.v 2 )v i + v. ri + vy 12 



dF 

dv\ 



O 



■ (4-7) 



We have replaced here the vectorial distance r[ by the un-primed notation ri, noticing that r[ 
is the dummy variable with respect to which the regularization proceeds (with this notation 
r' 2 is replaced by ri + y' 12 ). Following ( |3.35| ), the Lorentzian regularization [F] 1 is achieved by 
posing V = vi and v[ = 0, as well as y' 12 = y' 12 [yi2; Vi, v 2 ] and v' 2 = v 2 [yi 2 ; vi, v 2 ], where 
the latter functionals are defined in the general case by ( |3.29| ) and ( |3.30| ). It is convenient to 
obtain first an intermediate formula by setting V = vi and v[ = in ( |4.7| ), and by replacing 
into the terms that are already of order 1/c 2 the primed variables y' 12 and a' 1; a' 2 by the 
un-primed ones. Using also the identity ( |4.1|) , we arrive at 



[F] 1= [ F [n, r a + y' 12 ; v 1; Vl + v' 2 ] + _( Vl . ri )^^F 
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+ 



vi.ri, 



dF 
dy\ 



dF 



dy. 



dF 

dv{ 



dF 

'dvJ 



c- 



2 v l 



Vi.V 2 



^2 + (vi.yi 2 Ja 2 — - 



dF 

dvl 



1 / x r 1 i ti dF 

+ ^(vi.yi 2 )[-^l + ^ 



+ 



(4- 



where y' 12 and v 2 in the first term of the right side are given functions of yi 2 , Vi and v 2 
obtained by approximating (|3.29| ) and ( |3.30| ) to the first order; we find 



y n = yi2 + -o(vi.yi 2 ) 

c z 



-Vi + v 2 + 



v 2 



+ 



-0IV1.V2JV1 



bf\ + vi.v 2 



v 2 - (Va.yi2ja 2 



1 



(4.9a) 
(4.9b) 



(where the acceleration is equal to its Newtonian value). By inserting (|4.9|) into (|4.8|) and 
expanding to order 1/c 2 , it is easily seen that we cancel out exactly the two last terms in 
the right-hand side of (fOf), so that the result simplifies appreciably: 



F [n, r 2 ; vi, v 2 ] + — (vi.ri)^^F 



+ 



. l dy\ 2 dy\ 1 dv\ 



O 



Finally, we recognize in the right side the partial time-derivative, 

s dF 



d t F 

so that our final result writes 



dF , dF , dF 
dy\ 



F + ^(v x .W 



d t F + \v\diF 



2 dv\ ' 



(4.10) 



(4.11) 



(4.12) 



The result ( 4.12| ) displays the first relativistic corrections brought about by our 
Lorentzian regularization [F] 1 . As a check of the formula, let us apply it to the case of 
the special combination d{F — which, as we know from (|2.5|) , has no partie finie at the 
point 1 in the sense of the old regularization. This is no longer true in the sense of the new 
regularization. Using the equation (|4.12|) we find instead 



Ma 



v ni.vi 



AdtF 



O 



(4.13) 



The check consists of remarking that because of ( |2.5| ) we have (d[F' — 3^p L F') l = in the 

li 

rest frame of the particle 1, therefore the equation [d^F' — 3— rF'^ = must hold in any 
frame by definition of the new regularization. In the frame where the particle velocity is Vi 
we have r[ = n + ^(vi.r^vi + O(^) and d[ = + ^v\d t + ^v\v{dj + O(^). Inserting 
these relations into the previous equation, and using the fact that F is a scalar, we recover 
the formula ( |4.13|) after a short computation. 
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V. THE STRESS-ENERGY TENSOR OF POINT-PARTICLES 



With the Lorentzian regularization in hands, we make a proposal for the description of 
point-like particles in (post-Newtonian approximations of) general relativity. We recall first 
the general context of the problem. We want to solve the field equations of general relativity 
by means of analytic post-Newtonian series, with matter source describing appropriately 
defined point-particles. The stress-energy tensor of the matter source is supposed to be 
spatially isolated; we recall that, in this case, general relativity admits the Poincare group 
as a global symmetry. We assume the existence and unicity of a global harmonic coordinate 
system, defined by the gauge conditions 

<W = 0, (5.1a) 
hr = ^gsr-rr. (5-lb) 

where g^ u denotes the inverse of the covariant metric g^ u , and where g is the opposite of its 
determinant. The harmonic gauge conditions ( |5.1|) introduce a preferred Minkowskian struc- 
ture, with Minkowski metric given by r]^ = diag(— 1, 1, 1, 1) = r]^ v . Thus, the gravitational 
field can be described in harmonic coordinates by the Lorentzian tensor field propagat- 
ing on the Minkowskian background rf" . Similarly, one can think of the trajectories of the 
particles as accelerated world lines in Minkowski space-time. Subject to the conditions ( |5.1|) 
the Einstein field equations take the form of wave equations on the flat background, 

UhT = ^GgT^ + A" u [h,dh,d 2 h] , (5.2) 

where the flat d'Alembertian operator is given by □ = rf^d^dy. The right-hand side is made 
of the sum of the matter source term, with spatially compact support, plus the gravitational 
source term A Miy , given by a certain functional of the field variables h pa and its first and 
second space-time derivatives, and at least of second order in h. A consequence of the 
harmonicity conditions is that 

^(^ + ^)=0, (5.3) 

which is equivalent (through the contracted Bianchi identity) to the covariant conservation 
of the matter stress-energy tensor T^, 

W U T» U = , (5.4) 

the latter equation being in turn equivalent to 

du (V9 9x,T^) = \^gd x g^T^ . (5.5) 

In this section we regard the matter tensor T pv as a Lorentz tensor defined with respect to 
the Minkowski metric rj^ singled out by our choice of harmonic coordinates. 

To define a model for point-like particles, we follow essentially the derivation of the 
stress-energy tensor of test masses moving on a fixed smooth background (see e.g. Weinberg 
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30| p. 360). However, in the case of "self-gravitating" particles, we do not have a smooth 
background at our disposal, and the metric becomes singular at the location of the point- 
masses. Essentially, we shall propose the value of the (post-Newtonian) metric coefficients 
on each of the particles to be given by the Lorentzian regularization defined in Section [Til] . 
This entails supposing that the metric coefficients belong to the class of functions T . This 



is correct up to the 2PN order |[L9| ; however, we know that the expansion of the metric 



coefficients (in harmonic coordinates) near the particles, instead of being of the type ( |2.1|) - 
( |2.2| ), involve some logarithms of the distance to the singularities starting at 3PN order. 
It was shown that, at this order, the logarithms can be considered as some constants 
and included into the definition of the partie finie; moreover, they can be finally eliminated 
from the equations of motion by a change of coordinates. This suggests that we might 
consider more generally the logarithms as some constants, motivating our assumption that 
9fiu £ On the other hand, it is known ]2"2] , |2~3| ] that the constants Si and s 2 entering 



the partie-finie integral (2J3) must be adjusted in order that the equations of motion can 
be deducible from a Lagrangian, and in particular admit a conserved energy. For these 
reasons (presence of logarithms, equations of motion not directly admitting an energy), the 
following derivation of the stress-energy tensor for particles cannot be considered to be a 
rigorous proof. However, as we shall see, it is nicely consistent with the regularization, and 
its result satisfying. Our basic assumption is that the dynamics of the particles follows from 
the variation, with respect to the metric, of the action 



'particle Tfl\C 



+00 



dt ^-[^l^K + 1-2, (5.6) 



where = (c, dyi/dt) denotes the coordinate velocity of the particle 1 (we consider a two- 
body system, but the generalization to N bodies is immediate). The crucial point is that the 
value of g^y at 1 is assumed to be given by the Lorentzian regularization defined in Section 



(HJ. We vary the action ( |5.6| ) with respect to the metric, i.e. we imagine that g^ v G T 
is subject to an infinitesimal variation g^ v — > g^ v + Sg^ and compute the corresponding 
change in the action. However, we want the variation of the metric to correspond to the 
same matter system with two singularities 1 and 2. The evident and most natural way to 
ensure this is to suppose that Sg^ G T . Under the latter variation the regularized value 
of the metric at the point 1 undergoes the infinitesimal change [g^ u ]i — * [9nv\i + [&9nv\i- 
Therefore, the variation of the action (|5.6|) reads as 



Wparticle = \m X C / dt 11 = [Sg^ + 1^2. (5.7) 

) V _ i9fxr\l V l V l 



From the defining property ( |3.36| ) of the delta-pseudo-function PfA 1; we can re- write 
in the equivalent form 

r+oo v n v v 

^particle = \m x c / dt 1 1 = < Pf A 1; 5g^ > +1 «-> 2 . (5.8) 

J-oo V-Wi»l»' 
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Now, recall that the duality bracket is defined by the partie finie of the three-dimensional 
integral [cf (|2.8|)], so the latter expression can be cast into the standard form appropriate to 
the definition of a stress-energy tensor Tp a ^. ticle , namely 



/+oo 
dt < Vg T pMe » s 9n» > ■ ( 5 - 9 ) 
00 



The only difference with the standard definition is that the partie finie takes care of the 
divergencies at the positions of the particles. By comparing ( |5.8|) and (|5\9| ), we readily find 
that the corresponding stress-energy tensor density is given by 



V9 particle = rrnc f V \ = = PfAi + 1 <- 2 . (5.10) 
V [9pcrh v l v l 

The stress-energy tensor itself comes immediately from the rule of multiplication of pseudo- 
functions (|3.371) : 



This tensor takes the same form as the stress-energy tensor of test particles moving on a 
smooth background, except that the role of the background field is now played by the metric 
generated by the particles, regularized following the prescription ( |3.35| ) . Notice in particular 
that the factor 1/y/g inside the partie finie sign Pf should not be replaced by its reguralized 
value at 1 [see Q3.39Q 1. We propose the tensor ( |5.11| ) as a model of particles in the post- 
Newtonian approximation. From the product rules for pseudo-functions, we get the matter 
source term in the right-hand side of fl5.2|) as 



particle = rmc f lV ] = = Pf(V?AO + 1 - 2 . (5.12) 



The post-Newtonian iteration of the field equations in [p2] , p4]] is based on the latter expression 
of the matter source term. 

We now derive the equations of motion of the particle 1 from the covariant conservation 
of the stress-energy tensor ( p. lip : 



V,T p ^ rticle = 0. (5.13) 



Notice that thanks to the presence of the delta-pseudo-function, we know that the derivative 
is "ordinary" and satisfies the Leibniz rule in the sense of (|3.41|) . Thus, we can transform 



^ ^particle m ^ ne standard way and find that the equation ( |5.13j ) is equivalent, like in the 
case of continuous sources, to the alternative form 

du (V9 9x^T^ Iticle ) = \^gd x g^T^ n . cXe . (5.14) 



Then, we integrate ( |5.14| ) over a closed volume Vi surrounding the particle 1 exclusively. The 



role of the three-dimensional integral is played here by the duality bracket defined by 
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Let us denote by lvx the characteristic function of the volume Vi, such that lvi( x ) = 1 if 
x G Vi and lvi( x ) = otherwise [notably, ly 1 (y2) = 0]. Thus, we consider 

< d v {^fgg^ ^ tlcle ) , l Vl >=< \y/gd x g^ T p ^ rticlc , l Vl > . (5.15) 

(Though l-y-j does not belong to the class J 7 , it is locally integrable on M 3 and we know that 
the duality bracket applies on such functions as well; see 0.) The partial derivative d u in 
the left-hand side is split into a time-derivative and a space-derivative. Following the rule 
( P-40| ), the spatial derivative <9j is shifted to the right side of the bracket, where it applies on 
the characteristic function l Vl . Because of the presence of the delta-pseudo-function, the 
derivative of ly 1 is to be taken in an ordinary sense and is zero. Following the rule (9.7) in 
0, an analogous reasoning is valid for the time-derivative do = ^d t which can thus simply 
be put outside the bracket. Thus, we get 

^{< V9 9\»T^ rticle , l Vl >} =< ^Vgd x g^T^ ticle , l Vl > . (5.16) 

Next, we insert into ( |5.16| ) the specific expression ( |5.1(J| ) of the stress-energy density of 
particles. Because of the presence of the function ly 1 only the part corresponding to the 
particle 1 contributes, and we obtain 



d 
dt 



Finally, the effect of the brackets in both sides of the latter equation is to take the value 



at the point 1 in the sense of the Lorentzian regularization ( |3.35| ). Thereby our final result 
reads as 



The equations of motion of the particle 1 have the same formal structure as the geodesic 
equations of a test particle. In separate papers p2| , p4] , p5] ] we use (|5.18| ) to derive explicitely 



the equations of motion of the two particles at the 3PN approximation. 



ACKNOWLEDGMENTS 

We thank Misao Sasaki, Hideyuki Tagoshi and Takahiro Tanaka for stimulating discus- 
sions, and for the permission to reproduce their proof of the mathematical formula ( A29 ) 



at the end of the appendix. The support of the JSPS short-term program for research in 
Japan is gratefully acknowledged for a visit of one of us (L.B.) at the University of Osaka 
during which this work was begun. 



25 



APPENDIX A: SOLUTION OF THE EQUATION (|TT|) 



We are looking for the vector zi satisfying the equation 



zi = YiU - -»V.(x-zi' 



(Al) 



where yi(i) represents a given smooth (C°°) time-like trajectory and V a constant vector 
with norm |V| < c. Clearly, for a given trajectory, the solution z\ depends on the field point 
x as well as on time t. It was shown in the text after (|3 .15 ) that the application x — > z\ 
is contracting with fixed point yi. Here, let us look for the solution zi in the form of a 
function of the coordinates, 



Zi = zi(x,t) . 



(A2) 



From QAlp we compute the partial derivatives of z x with respect to t and x l , considered to 
be independent, and readily obtain 



dzi 
dx i 
dz\ 
~dt 



Vi-V 



dz\ 

dx i 



v a t--v.(x- Zl ; 



i + Iv*i 

c 2 dt 



vi i ~V.(x - z i' 



(A3a) 
(A3b) 



Contracting these equations with the vector V we can obtain the scalar products V.^ and 
V.%*, and use them back into ( |A3| ) with the result that 



dzi 
dx i 
dz\ 
~dt 



-iv. 



Vl 



V.vi 



Vl 



1 - 



V.vi ' 



(A4a) 
(A4b) 



where the velocity vi is evaluated at the instant t — 4-V.(x — zi). In particular, we find that 
zi must be a solution of the following first-order differential equation: 



dzi 

dx i 



c 2 1 dt 



(A5) 
) and 

(A6) 

necessarily satisfies the original equation ( |A1|) . Such a zi(x, t) being given, we perform in 
the equation (|A5|) the change of variables (x l ,t) — > (p\,Ti) defined by 



Conversely, let us prove that a vector Zl that (i) satisfies the differential equation 
(ii) admits yi(t) fixed point, i.e. is such that 

zi (yi (*),*) = yi(*) , 



pi = x i -z i 1 (*,t) 



Tl = t 



zi(x,t) 



(A7a) 
(A7b) 
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Using (|A5 ) it is easy to obtain the laws of transformation of the partial derivatives: 

d did 

+ -,V~ , (A8a) 



dp\ dx l c 2 dt 

where B % ^ denotes the matrix inverse of A\ = 5{ + jsVk-^- (i.e. A % -B\ = S l k ; in the case 
considered here where the velocities are strictly less than c the matrix A 1 • is a deformation 
of the unit matrix and thus admits an inverse). Now, under the change of variables (|A7| ) 
the differential equation ( |A5| ) becomes simply 

g-0, (A9) 

whose general solution is an arbitrary function of the time variable T\. Therefore, there must 
exists a trajectory Yi such that 

z 1 =Y 1 (r 1 )=Y 1 ^-lv.(x-z 1 )J . (A10) 

Imposing now that yi(t) is a fixed point for this solution zi in the sense of (|A6|) leads 
immediately to 

Yi(t)=yi(t), (All) 

so the equation (|A1|) is recovered exactly. Thus, solving (|A1|) is equivalent to solving the 
differential equation ( |A5| ) supplemented by the condition (|A6|). Notice that from ( |A1| ) or 
equivalently from ( |A5| )-(|A~6|) we find that Zi tends to the fixed point in the "non-relativistic" 
limit c — > +00, i.e. 

lim {z 1 (x,t)} = y 1 (t). (A12) 

c— >+oo 

This suggests to look for the solution Zi in the form of an infinite series of relativistic 
corrections of successive orders l/c 2n [from (|A5|) we know that z x is a function of 1/c 2 ]. 
Thus, taking also into account the limit (|A12|) , we pose 

+00 _ 

El n 
^Zi(x,t), (A13) 

n=l ° 



and we look for each one of the unknown coefficients Zi (x, t). By replacing the series (|A13|) 
into both sides of the equation (|A5|) and identifying the factors of the powers of 1/c 2 on each 
side we find, for any n > 1, 
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with the convention that Zi= yi(t). The equations ( |A14j ) are to be solved using the condition 
of fixed point yi \cf (|A6|)], which implies that, Vn > 1, 



The solution of (TAT§) - 
suppose that 



Zi ( yi (t),t) = 0. (A15) 
is found by induction over n. As an induction hypothesis 



n-l 

Zi = 



\n-l 



fn- 1)! V<% 



5 



n-2 



(A16) 



where ri = x — y 1; and where the partial time derivatives act on t keeping the space 
coordinate x fixed: for instance, dr\/dt = — Vi and dvi/dt = dvi/dt = a 1; where a x is the 
acceleration. Notice that (|A16 ) satisfies the condition (|A15|) because it involves n — 2 partial 
time derivatives while there is a factor (V.r!)™" 1 inside the brackets, so after differentiation 
there will remain at least one factor V.ri making the result be zero when x = y^. Inserting 
( |A16 ) into the right-hand side of (|A14|) we obtain the equation to be satisfied for the next- 
order coefficient, 



n 

dZi 



V, 



d 



n-l 



dx % 1 (n — 1)! \dt 
which can be re-written equivalently in the form 

n-l 



(V-r^-Vi 



(A17) 



dZi 

dx l 



d_ ((- 

dx 



dtj 



(V. ri )"vi 



(A18) 



showing that the most general solution is necessarily of the type 

(V.nyVxl + C(t) , 



n 

Zi= 



nl 



- 



dtj 



(A19) 



where C(t) denotes an arbitrary vector depending only on time t. However, this vector must 
be zero on account of the fact that the result should be zero when x = yi. Therefore we 
have proved by induction that 



nl 



-) 

dtj 



n-l 



so the vector zi solving at once (|A5| ) and ("KB), or equivalently 
rather interesting infinite series 



zi = yi 



+oo 
n=l 



d_ 

c 2n n\ V dt 



n-l 



;v. ri )" Vl 



(A20) 

), takes the form of the 
(A21) 
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which constitutes the solution needed for our work in Section [II. Furthermore, subtracting 
x from this solution and contracting with V we obtain after a short calculation the quantity 
Ti which was defined in ( |A7b| ): 



n 



<+E 



_N3 (— ) n ( d N 



n=l 



c 2n n\ \dt 



(A22) 



Now, recall that the latter quantity T\ is such that zi = yi(ri). Therefore, we see that 
we can find an alternative expression of the vector Zi by inserting into yi(ri) the series 
expansion ( |A22j ) found for t\. Using an infinite Taylor expansion we are led to 



zi = yi 



+00 

+E 

p=0 



(p + 1)! dtp 



+00 

E 

, n=l 



c 2n n! \<9t 



5 



n-1 



p+1 



:v. ri ; 



(A23) 



Each of the terms is composed of p + 1 sums; accordingly we introduce p + 1 summation 
indices n±, . . ., n p , n p+ i so that 



z-i 



^oo +00 



E-E E 



(itP c 2(m+...+n p+ i) 
p=0 ni=l n p =ln p +i=l 



X 





ni— 1 


"(V.r!) ni " 






nil 


pose 


jfe = 


m + . . . - 



n,- lr 



"(V.ri) 71 *^ 




"(V.ri) n f+ 1_ 






n p+ i! 



(A24) 



commutations of summations to arrive at 

-- 'X 

C 



zi = yi 



d p v l 



E-E 



fe=i 



(p+1)! g^p 

p=0 v ^ 7 m=l n p =l 



X 



in which n 





^11— 1 


"(V.ri) ni " 






Til! 



\dtj 



n„-lr 



"(V-ri)""" 


/ 9 yp+i-i 


"(V.ri) n f+ 1_ 


n p ! 




n p+i ! 



P+1 



p+i, 
t=j 



71; 



(A25) 



1) (with 1 < j < p). We must identify 



the latter complicated expression with the simpler form of the vector z x given by ( |A21[ ). 
From identifying the powers of 1/c 2 in both expressions we immediately obtain 



fc-i 



(V-ri) 

jfe! 



-Vi 



fc-l 1 ,„ qi 1v 

El crvi \ ^ \ + 

(p + i)! dtp 2s ' - 2s 

p=0 K ^ ' m=l n»=] 



yi-i 


"(V.ri)" 1 " 


/ q yp-i 


"(V.ri) n f" 


/ 3 \ "p+i-i 


"(V.ri) n f+ 1 " 




ni! 


"UJ 


n p \ 


UJ 


n p+ i! 



(A26) 



Finally, from using the binomial formula for the derivative of a product, we can identify 
in each side of the latter equation the coefficients of each d p \i/dt p , and we arrive at the 
relation, valid for any p and any k > p + 1, 



29 



91 



• • • ( m 



?ll=l 



n„=l 





"(V. ri ) ni " 


/ 9 y P -i 


"(V.ri) n " _ 


/ 9 yp+i-i 


"(v.n)™^ 1 " 


w 


Till 


"' UJ 




UJ 





(p+l)(fc-l)! 



(A27) 



The latter relation actually represents a quite general mathematical formula because we 
have specified nothing about the scalar product V.ri- Therefore, the relation ( |A27| ) holds 
in fact in the case of an arbitrary sufficiently differentiable function f(t), so 



E-E 



rti=l 



n„=l 



dt 



nil 



n K -l 



f n p 



n 



d y^ 1 - 1 

dt) 



f 



n p+ i\_ 



(p + l)(fc- 1)! /_d 
(Jfe-l-p)! Vdt 



fc— p— i 



k\ 



(A28) 



The equivalence obtained above between the formula ( [AT] ) and the differential equation 
together with the auxiliary condition (|A6|) shows indirectly that the mathematical 
formula (|A28[) is correct. However, a direct proof of this formula has been found by Tanaka, 
Sasaki and Tagoshi. Here we reproduce their proof in the particular case where p — 1, so 
that qi — k — 1 and n 2 = k — n (where n = rij), in which case the formula reads, for any 
k > 2, 



jfc-i 

E 

n=l 



n-1 



nl 



k—n—lr 



k—n 



(k-n)\ 



2(k 



k - 2 r jk-[ 



k\ 



(A29) 



We replace f(t) in ( |A29|) by its Fourier transform, f(t) = f*™ ^e lu}t f(u), and readily find 
that in order to prove the formula (|A29 ) it suffices to prove the statement that the equation 



k-l 

E 

n=l 



LO(\ + UJ 2 + . . . + UJ r 



n-1 



LO n+ i + . . . +Uk) 



k—n—1 



2(k - 1) (uJi + u 2 + ■ ■ ■ + uj k 



k-2 



(A30) 



holds identically for any family of real frequencies u>i, u> 2 , ■ ■ ., 0Jh- Most importantly, the 
parenthesis around indices in the left side of (|A30| ) indicate the complete symmetrization 
over the k frequencies Ui, . . ., uj^ [in addition, ( ) denotes the binomial coefficient]. Let us 
single out one of the frequencies, for instance u>k, and re- write (|A30|) in a form involving an 
explicit symmetrization over the other k — 1 frequencies, U\, . . ., uJk-i, only: 



fe-i 

E 

n=l 



k-l 

n 



W(i + ■ ■ ■ 



n-1 



^>n+l + • • • +UJk-l) + ^k 



k—n— I 



(k — 1) [ut +u 2 + . . . + u k 



k-2 



(A31) 
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(in which we have simplified a factor 2 in both sides of the equation). Furthermore, let us 
replace in the latter formula ujk by some sum Uk + . . . + Uk+ S , and symmetrize over the whole 
set of frequencies u>i, . . ., oj^+s- This yields, for any s, 

E' 1 /jfe-l\ / \ n ~V \fe-n-l 

I J \ U(i + . . . + U n J UJn+1 + ■■■ +U k +s) J 

n=l ^ ' 

= (A;-l)(wi + w 2 + ... + w fc+a )* 2 . (A32) 
Now we prove that the equation ( A30 ), or equivalently ( A3 1| ) , is true by induction on the 



integer k. Therefore, our induction hypothesis is that (|A31|) is correct for any k < K, and 



from this we want to show that it is correct again for k = K + 1. Note that from our 
induction hypothesis we know that ( |A32|) is also correct for any k < K and any s. Consider 
the sum defined by the left side of (|A31|) in the case where k — K + 1, say 

Sk+x = ( j + • • • + uj n j (uj n+1 + . . . + u K) + uJ K +i) , (A33) 

n=l V n / 

where we recall that one of the frequencies, i.e. uk+i, is "artificially" singled out. However, 
Sk+i is also given by half the left-hand side of ( |A30| ) and is symmetric in ux, . . ., lor+i- 



We want to show that Sk+i is equal to the right-hand side of ( |A31| ) with k = K + 1. To 
this end, we transform Sk+i with the help of the binomial formula, and obtain after a short 
calculation 

S K+1 = 2^ n (K^Bl l^y n )[U(1 + ---+Un) [u n + 1 + ...+U J K )j 

n=l 



1=0 



(A34) 



Now we have two sums over I and n, and it is easy to recognize that the second sum, over 
n, can be simplified as soon as I > 1 by means of ( [A32j ) which is correct by induction under 
the condition that k < K and for any s. Posing K — I = k — 1 and k + s = K we see that 
this condition is realized if and only if I > 1. After simplification we find 



( \K-\ 

S K +i = K\u 1 + ... + uk+i) + ^k+i (^i, • • • , uk) , (A35) 

where the first term is the result we want to obtain, and where the second term is a certain 
function of the frequencies u\, . . ., ujk but which does not depend on luk+i- The expression 
of ^ k+i is given for completeness as 

K 



^K+l = ^2 \ {^i 1 + ' ' ' + w n) (un+1 + • • • + UK)) ~ K (u X + . . . + UJ K ^j 

n=l V ^ / 

(A36) 
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Now we use the fact that Sk+i is actually fully symmetric with respect to the K + 1 
frequencies ui, . . ., uk+i- Therefore the function ^k+i must be a pure constant, independent 
on any uo n . Furthermore, we know also that Sk+i is a homogeneous polynomial of degree 
K — 1 in all the u\, . . ., lok+Ii so this constant must in fact be zero: ^^+1 = 0. Finally we 
are able to conclude on the desired result, 



Incidentally, notice that the equality ^k+i — is itself a consequence of the same mathe- 




(A37) 



matical formula, since it follows from setting k = K + 1 and posing ujk+i 



in flMTp . 
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